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(N 

^H We present several new spatially localized equilibrium and traveling-wave solutions of 

^-( plane Couette and channel flows. The solutions exhibit strikingly concentrated regions of 

vorticity that are flanked on either side by high-speed streaks. For several traveling-wave 
solutions of channel flow, the concentrated vortex structures are confined to the near- 
P\| wall region and form particularly isolated and elemental versions of coherent structures 

in the near-wall region of shear flows. The solutions are constructed by a variety of meth- 
ods: application of windowing functions to previously known spatially periodic solutions, 
^L^ continuation from plane Couette to channel flow conditions, and from initial guesses 

''O obtained from turbulent simulation data. We show how the symmetries of localized solu- 

r^ tions derive from the symmetries of their periodic counterparts, analyze the exponential 

y^ decay of their tails, examine the scale separation and scaling of their streamwise Fourier 

modes, and show that they develop critical layers for large Reynolds numbers. 

O 

^ 1. Introduction 

r^ Over the last twenty years the computation of invariant solutions of the Navier-Stokes 

Oh equations, or "exact coherent structures," has opened a new approach to understanding 

the dynamics of moderate- Reynolds unsteady flows, an approach which promises to pro- 

,_^ vide a long-hoped-for bridge between dynamical systems theory and turbulence. Unlike 



< 



a 






^ previous derivations of low-order dynamical models of unsteady flows (Lorenz (1963); 

C^ Aubry et al. (1988); Holmes et al. ( 1996[ )), the invariant-solutions approach forgoes low- 

cn 



d projections and simplified models and instead takes a fully-resolved direct numerical 
simulation as a quantitatively accurate finite-dimensional approximation of the Navier- 
^^ Stokes equations. The fully-resolved simulation is then treated as very high-dimensional 

^^ dynamical system. The first step in analysis of a dynamical system is the computation of 

«y~, its invariant solutions: its equilibria, homo- and heteroclinic orbits, and periodic orbits. 

,-H In the dynamical systems view of turbulence, the equilibrium solutions correspond to 

ILjl steady states of the fluid flow, periodic orbits correspond to states of the fluid velocity 

field repeat themselves exactly after a finite time, and homo- and heteroclinic orbits 
correspond to dynamic transitions between equilibria or periodic orbits. For flows with 
homogeneous spatial directions, such as pipes and channels, continuous symmetries in 
the equations of motion allow relative invariant solutions, e.g. traveling waves (relative 
equilibria) and relative periodic orbits. 

The simplest invariant solutions of fluids are the classical, closed-form steady states of 
the Navier-Stokes equations, for example, the parabolic laminar flow profile of pressure- 
driven channel and pipe flow, or the linear laminar solution of plane Couette flow. So- 
lutions such as these have special cancellations which make it possible to represent the 
exact solution of the nonlinear system with a finite set of simple functions. For example, 
for the laminar solution of channel flow, the nonlinear term vanishes and the solution can 
be represented exactly as a 2nd order polynomial in the wall-normal variable. However, 
if we consider the problem from the perspective of faithful, very high-dimensional finite 
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discretizations, invariant solutions are the solutions of a nonlinear algebraic or differen- 
tial equations in 0{\Q'^) to O(IO^) free variables. Only a few very special solutions (the 
classical ones) will involve few enough modes to be expressible in closed form, and most 
will involve nonlinear coupling between large numbers of nonzero variables. Compared 
to the classical closed-form solutions, these computed invariant solutions are typically 
unstable, fully three-dimensional, fully nonlinear, distant from the smooth laminar flow 
solutions, and involve most if not all the available modes of the numerical representation. 
Determination and specification of such solutions is necessarily numerical. 

The practical feasibility of finding such high-dimensional nonlinear solutions of Navier- 



Stokes was first demonstrated by Nagata ( 1990 1, who computed an unstable 3-dimensional 



nonlinear equilibrium solution of plane Couette flow at a Reynolds number above the on- 
set of turbulence, using a 589-dimensional discretization. The same equilibrium solution 
was found independently and analyzed in greater precision and detail by [Clever fc Busse] 
([l992j and Waleffe (1998( 2003). A large number of equilibria and traveling waves of 



plane Couette and pipe flow have since been found (Nagata (1997); Schmiegel 



Faisst & Eckhardt (20031; Gibson et al. (2009)), a few of channel flow (Waleffe 



Itano & Toh (2001)), and in other flows such square duct flow (Okino et al. 



(1999); 



(2001); 



(2010), 



Wedin & Kerswell (2009)). Periodic orbits have been calculated for plane Couette flow 
(I Kawahara fc Kida| (|2001[ ); |Viswanath| ( |2007[ ); [Cvitanovic &: Gibson] (|2010[)), p i pe flo w 
lOuguet et aL| ( |2008| ), and 2D Kolmogorov turbulence ( [Chandler fc Kerswell] (|2013[)) , 
and hetero- and homoclinic connections for plane Couette flow (Gibson et al. (2008); 



Halcrow et al. (2009); van Veen fc Kawahara (2011)). Improved numerical methods and 
more powerful computers now allow the computation of solutions with as many as 10^ 
free variables. High-resolution calculations have shown that discretization errors con- 
verge toward zero as resolution is increased, demonstrating that the numerical solutions 
are precise approximations of true solutions of the Navier-Stokes equations, rather than 
artifacts of discretization. High-resolution calculations have also allowed accurate com- 
putation of solutions with flne spatial structure, such as periodic orbits that exhibit 



turbulent "bursting" phases (Viswanath (20071; Cvitanovic fc Gibson (2010)). 

Just as in low-dimensional dynamical systems theory, the importance of these invari- 
ant solutions stems from the organization they impose on the state space dynamics. In 
particular, dynamics in the neighborhood of (relative) equilibria and periodic orbits is 
governed to leading order by the linearization about these solutions, and the eigenvalues 
of the linearized dynamics reveal the local character of the state space flow and the di- 
mensionality of each solution's unstable manifold. For shear flows at moderate Reynolds 
numbers and in closed or small periodic domains, most known invariant solutions have 
a positive but remarkably small number of unstable eigenvalues, and correspondingly 
low-dimensional unstable manifolds. For example, the equilibrium solution of plane Cou- 
ette flow developed by Nagata, Busse, Clever, and Waleffe (hereafter termed the NBCW 
equilibrium) has a single unstable eigenvalue (Wang et al. ( 2007| ), and the periodic orbit 
solutions of Viswanath (20071 have between 1 and 11 unstable eigenvalues. This low di- 



mensionality of instability is a crucially important result. It suggests, as long suspected, 
that moderate-Reynolds flows are inherently low-dimensional -at least for small conflned 
domains. It further suggests that the temporal dynamics of such flows results from a 
relatively low-dimensional, chaotic but deterministic walk between the flow's unstable 



invariant solutions, along the low-dimensional network of their unstable manifolds ( Gib 



son et al. (2008)). Moreover, the coherent structures often observed in such flows can be 



understood as resulting from close passes to these unstable invariant solutions, at which 
the Navier-Stokes equations balance exactly -an idea well-expressed by Waleffe's term 



"exact coherent structures" (Waleffe (2001)). Indeed, a key feature of the NBCW invari- 
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ant solution is that it captures structure commonly observed in shear flows in the form 
of wavy rolls that support alternating streaks of high and low streamwise velocity (see 



2.2 for further discussion). We refer the reader to the Kawahara et al. (2012) review 



article for an excellent overview of research in this area. 

On the other hand, a significant weakness of the invariant-solutions approach to date 
is the assumption of idealized computational domains, for example, small periodic "min- 



imal fiow units" (Hamilton et al. (1995)), just large enough to contain a single coherent 
structure. For example, most of the work analyzing dynamics via invariant solutions has 
been done in plane Couette flow in doubly-periodic boxes with stream- and spanwise 
aspect ratios around 2 or 3. Small periodic or confined domains are an understandable 
simplifying assumption for the demanding computational problem of computing invari- 
ant solutions. However, it is also valid to criticize such assumptions as unphysical and 
limiting the research to explorations of temporal complexity of periodic spatial structures, 
as opposed to the full spatio-temporal complexity of turbulence in extended domains. The 
assumption of periodicity also complicates coordination of theory and experiment. Close 
passes to unstable traveling waves with axial periodicity have been have been detected in 



experimental pipe flows (Hof et al. (2004); de Lozar et al. (2012)), but the effort to match 
experiment and theory would be greatly aided if the assumption of periodicity in the com- 
putations could be relaxed. Lastly, the most prominent example of coherent structures 
in fluids, and the initial motivating problem for dynamical-systems approaches, are the 
lambda and hairpin vortices that form in the high-shear region near the walls of channel 
and boundary layer flows. These appear in coupled systems of coherent structures that 
are individually localized in the wall-normal direction as well as the two homogeneous 



directions (Adrian (2007)) 



The first step in addressing the weakness of small, idealized domains is to compute 
localized invariant solutions on extended domains. Several papers have made promis- 



ing steps in this direction. Schneider et al. (2009) computed the first spatially localized 
equilibrium and traveling waves of Navier-Stokes, by an "edge-tracking" algorithm for 
plane Couette flow in a streamwise-periodic but spanwise-extended domain, and showed 
that these solutions were spanwise-localized forms of the spatially periodic NBCW solu- 
tion, which exhibit exponential decay towards laminar flow in the spanwise coordinate. 
Another computation in the same paper produced an unsteady, time-varying state with 
exponential localization in both span- and streamwise directions. However, this state ap- 
pears to wander chaotically and so is not an invariant solution. A similar time-varying 
doubly- localized state was found by Duguet et al. (20091. Schneider et al. ( 2010|) demon- 
strated a number of interesting connections between the localized solutions of Schneider 



et al. ( 2009 ) and localized solutions of the Swift-Hohenberg equation. Both systems ex- 



hibit homoclinic snaking, a process by which localized solutions grow additional structure 
at their fronts via a sequence of saddle-node bifurcations in a continuation parameter 
(see §4.2 for further discussion). This is an intriguing connection, as Swift-Hohenberg 



is a key model equation in the theory of pattern formation (Hoyle (20061), for which 
localization is comparatively well-understood (Burke & Knobloch (2007)). 

The broad purposes of the present paper are (1) to further extend the invariant- 
solutions approach to turbulence to spatially extended flows, towards the long-term goals 
of addressing spatio-temporal complexity in turbulence and aiding the effort to verify in- 
variant solutions in experiment, and (2) to begin an effort to capture spatially isolated 
coherent structures that occur near the walls of shear flows. The specific results and 
organization of the paper are as follows. In §[2] we show that the spanwise localized so- 
lutions of Schneider et al. (2010 2009) are not anomalous, but that localized versions 



of other spatially periodic solutions exist and can be constructed easily by a windowing 
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and refinement method that, unUke edge-tracking, puts no restrictions on the number 
of the solution's unstable eigenmodes. We show how the symmetries of localized solu- 
tions result from the symmetries and phase of the underlying periodic solution. In §[31 to 
further develop the invariant-solutions approach in experimentally accessible flow condi- 
tions, we construct localized traveling wave solutions of channel flow, by windowing and 
refining periodic solutions obtained by continuation from plane Couette conditions and 
by searching among turbulent simulation data. In doing so we find particularly intriguing 
traveling-wave solutions of channel flow whose vorticity is concentrated in the near-wall 
region, in spanwise and wall-normal localized structures that resemble lambda vortices in 
streamwise-developing channel flows. In §|4]we analyze the tails of the localized solutions 
and show that they decay exponentially to laminar flow at the rate determined solely by 
the streamwise wavenumber of the solution, with far-field structure that is independent of 
the details of the core region. We examine scale separation and scaling in the streamwise 
Fourier harmonics and development of critical layers at large Reynolds numbers. 



2. Equilibrium solutions of plane Couette flow 

2.1. Equations of motion 

Plane Couette fiow consists of an incompressible fiuid confined between two parallel rigid 
plates moving in-plane at a constant relative velocity. The x = {x, y, z) coordinates are 
aligned with the streamwise, wall-normal, and spanwise directions, where streamwise 
is defined as the direction of relative wall motion. We assume a computational flow 
domain fl — [~Lx/2, Lx/2] x [—h, h] x [—Lz/2, Lz/2] with periodic boundary conditions 
in X and z and no-slip conditions at the walls y = ±/i. We restrict our attention to 
streamwise-periodic velocity fields and L^ chosen to match the streamwise wavelength. 
In the spanwise direction, we choose L^ either to match the spanwise wavelength of a 
spanwise periodic field, or to a large value that approximates a spanwise-infinite domain. 
We decompose the total velocity and pressure fields into sums of a laminar base flow 
and a deviation from laminar: Utot(x, t) = u(x, t) + U{y) e^ and ptot — p(x, t)+ x dP/dx, 
where dP/dx is a fixed constant specifying an imposed mean pressure gradient. For plane 
Couette flow we will consider only the case dP/dx = 0, for which the laminar solution is 
U{y) — Uy/h, where U is half the relative wall speed. After nondimensionalization by 
U, h, and the kinematic viscosity z/, the Navier-Stokes equations for plane Couette flow 
can be written 

-^ + U^+vU' e^ + u-yu = -Vp-\ V^u, V-u = (2.1) 

ot ox Re 

where Re = Uh/v and the velocity components are u(x, t) = [w, u, u;](a;,2/, z, t). In this 
decomposition the plane Couette laminar solution is U{y) — y, dP/dx = 0, u = 0, and 
p = 0. From here on we refer to u as velocity and Utot as total velocity, and we note that 
u has Dirichlet boundary conditions at the walls. 



The Navier-Stokes equations (2.1 1 with plane Couette conditions and y-Dirichlet, x, z- 
periodic boundary conditions are invariant under any combination of rotation by tt about 
the z axis, reflection about the z = plane, and finite translations in the x and z 
directions. The generators of the plane Couette symmetry group are thus 



t{Ax,Az) 



[u,v,w]{x,y,z) -)■ [-u,-v,w]{-x,-y,z) 

[u, V, w] {x, y, z) -^ [u, V, -w] {x, y,-z) (2.2) 

[u, V, w]{x, y, z) — )■ [u, w, w\{x + Ax, y.,z + Az). 
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Figure 1. Three spatially periodic equilibria of plane Couette flow, (a) NBCW lower 
branch, (b) EQ7 lower branch, and (c) EQ8 (upper branch of EQ7), all at Re — 400 and 
Q,7 = 1,2, where a and 7 are the streamwise and spanwise wavenumbers. The visualizations 
show 3D isosurfaces of signed swirling strength at s = ±0.09 in green/blue (see text). Isosur- 
faces of streamwise velocity indicating high-speed streaks are shown in red at (a) u = 0.3, (b) 
u = 0.15, and (c) u — 0.2. The back plane shows contours of streamwise velocity at levels (a,c) 
±0.03 -[1:2:15] and (b) ±0.03 -[1:2:5] in Matlab notation with dashed/solid lines indicating 
negative/positive values. Note that by symmetry both solutions have streaks of equal magnitude 
and opposite streamwise velocity; these are not shown in the isosurfaces to reduce clutter, but 
they are indicated by the negative/positive symmetry of the backplane contours of streamwise 
velocity. The origin is at the center of the box. 



We use the standard group-theory notation (. . .) to indicate the group generated by a set 
of group elements; thus the symmetry group of plane Couette flow is {axy, (Jz, t{^x, Ax)). 
For each subgroup of this group, there is a subspace of velocity fields that is invariant 
under the equations of motion. That is, if a velocity field u(x, 0) satisfies u = cru for 
each symmetry cr in a given subgroup, u(x, t) will satisfy the same symmetries for all 
time. Invariant solutions of the equations of motion naturally lie in these subspaces. For 
example, with appropriate choice of the origin, equilibrium solutions of plane Couette 
typically satisfy [m, u, w](a;, y, z) = [—u,—v,—w\{—x,~y,~z) or related symmetries in- 
volving Gxy and Gz, since these symmetries requires the velocity field to vanish at the 
origin, and so prevent drifting in x and z. 

2.2. Spatially periodic solutions: EQl/NBCW, EQ7/HSV, and EQ8 
Figure IT] shows a visualization of three spatially periodic equilibrium solutions of plane 



Couette flow: the well-known "lower branch" solution of Nagata (1990), Clever & Busse 
( [l992| , and |Waleffe| ( p98] ) (N BCW, called EQl inpbson et a/.| ( |2"009f )), and the "Hair- 
pin Vortex Solution (HVS)" of Itano & Generalis (2009) and discovered independently as 



EQ7 in Gibson et al. (2009). Henceforth we refer to these as NBCW and EQ7. EQ8 is the 



upper branch of the the HSV/EQ7 solution. The NBCW solution is well-known not only 
as the flrst known exact nonlinear solution to the Navier-Stokes equations, but also for a 
number of remarkable characteristics, which we outline briefly here. The NBCW solution 
captures precisely, in the context of plane Couette flow, the roll-streak structure that 
seems ubiquitous in shear flows ranging from Taylor-Couette to the turbulent boundary 
layer, and consequently forms an example of an exact instantaneous balance between the 



three cycles of Waleffe's self-sustaining process for shear flows (Waleffe (1997 1998)) 



The NBCW solution has also served as a starting point for recent efforts to formulate a 
dynamical-systems theory of turbulence. At Reynolds numbers above the onset of tur- 
bulence, the NBCW solution lies between the laminar solution and the chaotic turbulent 



region of state space. It has a single unstable eigenvalue (Wang et al. (20071, so that its 



stable manifold forms a boundary between states that decay to laminar flow and states 
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that grow to turbulence (iSchneider et al. (2008)). The comparatively low viscous shear 
rate of the NBCW solution suggests that it might be feasible to implement a control 
strategy to stabilize this single unstable direction and obtain savings in the wall driving 



force compared to the turbulent flow (Wang et al. (2007)). The NBCW solution also 



has been shown to have well-defined asymptotic structure in the limit of large Reynolds 
number, which can be exploited to form a reduced system that accurately captures the 



structure of the solution over a wide range of Reynolds numbers (Wang et al. (2007); 
pi all fc S her win (20101). The asymptotic structure and reduced system are particularly 
relevant for this work, since it seems likely that any analytic understanding of localization 
in Navier-Stokes will be more easily developed in the context of a reduced system. 
The EQ7 solution has been conjectured to be related to hairpin vortices frequently 



observed in the turbulent boundary layer (Itano & Generalis (2009)). The hairpin shapes 
in the visualizations of Itano & Generalis (2009) are certainly suggestive; however, some 



caution should be in order, as these visualizations of vortex lines cannot be compared 
directly to visualizations that highlight the magnitude of vorticity, such as Q criterion, 
lambda criterion, and swirling strength. Figure [l] shows NBCW, EQ7, and EQ8 visual- 
ized with signed swirling strength isosurfaces to show the roll structure and streamwise 
velocity isosurfaces to show high-speed streaks. The swirling strength at x is defined as 
the magnitude of the complex part of the eigenvalue of the velocity gradient tensor Vu(x) 
(Zhou et al. (1999)). We chose swirling strength over other measures of fluid circulation 
such as Q criterion because it most clearly identified in 3D isosurfaces the regions of 
highly concentrated circulation that are apparent in 2D quiver plots such figure |4] Since 
the invariant solutions in this paper have elongated regions of concentrated circulation 
nearly aligned with the x axis, we attached a ± sign to the swirling strength that indicates 
clockwise/counterclockwise circulation with respect to the positive x axis. 

2.3. Construction of localized initial guesses by windowing 

The localized equilibria and traveling waves of plane Couette flow described inl Schneider| 
et al. (20091 and Schneider et al. (2010) (hereafter SGBIO) are spanwise-localized ver- 



sions of the spatially periodic NBCW solution. These localized solutions are comprised of 
a core region that closely resembles the periodic NBCW solution, weak tails that decay 
exponentially towards laminar flow, and a transitional region between the core and tails. 
This form suggests that new localized solutions might be found by imposing a similar 
core-transition-tail structure on other known spatially periodic solutions, and then re- 
fining these initial guesses with a trust-region Newton-Krylov solver. The rough form of 
this desired structure can be imposed on initial guesses by multiplying a known spatially 
periodic solution, expressed as a perturbation on laminar flow, by an even positive win- 
dowing function W{z) that is nearly unity over a core region \z\ < a, decreases smoothly 
and monotonically to nearly zero over a transition region a ^ \z\ < a + b, and vanishes 
as \z\ — )■ oo, followed by projecting the resulting field W^(z)u(x) onto the divergence-free 
subspace. We have found that with a robust trust-region Newton-Krylov solver, the pre- 
cise details of the windowing function and the projection are unimportant, and that the 
only important details are smoothness and the widths of the core and transition regions. 
One choice that suffices is the windowing function 



Wiz) ^ \ 



1 + tanh 



6(a 



l+tanh(^(^+3 



(2.3) 



This W{z) behaves as desired: it is even, smooth, monotonic in \z\, satisfies 0.995 < 
W{z) < 1 for \z\ < a and < W{z) < 0.005 for \z\ > a + b, and it approaches 
zero exponentially as \z\ — ^ cx). W{z) is specified in this particular form because we 
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found that the most important factor in producing a good initial guess was the size 
and location of the transition region, which are specified by the parameters a and b. A 
sufficient projection is to apply W{z) to the streamwise and spanwise components of 
velocity and reconstruct the wall-normal from the divergence-free condition. That is, let 
u = [u,v,w] be a ^-periodic solution expressed as a perturbation over laminar flow. An 
initial guess for a z- localized solution Ug = [ug,Vg,Wg] can be constructed by setting 
Ug — Wu, Wg — Ww, and reconstructing Vg from V • Ug = and boundary conditions. 

It is worth emphasizing that this localization procedure is rather crude. By construc- 
tion, the initial guess should nearly satisfy the Navier-Stokes equations in the core region 
and the tails -nearly but not exactly because the guess merely approaches the laminar 
solution for large z, and because the nonlocal effects of the divergence-free projection 
and pressure will corrupt the balance of terms that one would otherwise expect in the 
core region where W{z) is very nearly unity. In the transition region, however, there is 
no reason to expect that the velocity field that smoothly interpolates between tails and 
core will be close to satisfying Navier-Stokes. The quality of these initial guesses, thus, 
depends entirely on the robustness of the solver used to refine the initial guess into a 
solution. In particular, within the transition region the initial guess is too far from sat- 
isfying Navier-Stokes to be refined to an exact solution with a straight Newton method. 
A so-called globally convergent search method is instead required. We use a "hookstep" 
trust-region modification of the Newton search method coupled with a Krylov-subspace 



method (GMRES) for solution of the Newton-step equations, following Viswanath (2007 



2009). 



2.4. Localization and symmetry 

The symmetries of a desired solution are important both in determining solution type 
(e.g. equilibrium versus traveling wave) and for reducing the search space, which im- 
proves the speed and robustness of the search. The appropriate symmetries for spanwise- 
localized solutions are determined as follows. We begin with a spanwise periodic solution 
with a known set of symmetries. Multiplying that solution by a nonperiodic windowing 
function W{z) breaks any of these symmetries that involve z periodicity. Symmetries 
that do not involve z periodicity are preserved through the localization procedure and 
form the symmetry group of the localized guess. Note further the symmetry group G of 
a periodic solution u transforms by conjugation to tGt^^ when the solution is phase- 
shifted to TU, and that the localization procedure will break and preserve the different 
symmetry groups of u and ru differently. Thus it is possible to construct localized guesses 
with different symmetry groups by applying the windowing function to the same periodic 
solution in different spatial phases. 

To illustrate, we show how the symmetries of the equilibrium, traveling wave, and rung 



solutions of PCF in Schneider et al. { 2010 ) arise from localizing the NBCW solution in dif- 
ferent spatial phases. For compactness in what follows let Tx — t{£x/2, 0), r^ — t(0, iz/2), 
and Txz = T^Tz- In the spatial phase of Waleffe ( |2003 ), the (£2,,i!z)-periodic NBCW so- 



lutions have symm etry group {Txaz,Tza xyz) = {e,Txaz,Txz(Txy,Tzaxyz}, which is the S 
symmetry group of Gibson et al. { 2009 1 (hereafter GHC09) . If] The localizing procedure 



above sets [ug, Wg]{x, y, z) = W{z)[u^ w\{x^ y, z) and determines Vg from incompressibility. 
A simple series of substitutions shows that the first symmetry is preserved under local- 
ization, Ug — TajiTzUg, but the second and third symmetries are not: Ug ^ T^z'^xy'^g 
and Ug 7^ TzUxyzVig. Intuitively, since the windowing function W{z) is constant in 
X and y and even about z — Q but not periodic in z, windowing preserves the z- 



t Note that in Gibson et al. (20091, the y subscript on a^y was suppressed. 
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Figure 2. Construction of localized initial guesses by windowing. (a,b) The spa- 
tially periodic EQ7 solution of plane Couette flow, in two spatial phases. Three copies of the 
a, 7 = 1,2 periodic solution at Re = 400 are shown in the z € [— 37r/2, 37r/2] subset of the full 
Lx,Lz = 2tt, Gtt computational domain, with the z = plane shown bisecting the box to high- 
light the z symmetries. The solutions in (a) and (b) are related by a quarter-wavelength phase 
shift in z. (c,d) Initial guesses for localized solutions with different symmetry groups, produced 
by multiplying (a,b) by windowing functions W{z) (heavy Hne in front planes). Isosurfaces of 
signed swirling strength at s = ±0.12 are shown in (green/blue) and contours of streamwise 
velocity in the back plane at contour levels u — ±0.03 • [1:2:9]. 



reflection, x-translation symmetry r^az of the NBCW solution, but not its T^z^xy or 
Tzf^xyz symmetries, which both involve z periodicity. The sole preserved symmetry, 
T^Ozi is in fact the symmetry of the localized traveling wave reported in SGBIO, i.e. 
[wg, tig, Wg](a;, J/, z) — [ug, Wg, —W^(x + ^x/2, y, —z). Refinement of this initial guess by a 
search method that respects symmetry results in a traveling-wave solution with the same 
symmetries. 

The same localization process on a shifted NBWC solution produces an initial guess 
with the symmetry of the localized equilibrium solutions of GHC09. Shifting the NBCW 
solution by a quarter-wavelength in 2, i.e. r = r(0,^z/4), thus changes its symmetry 
group by conjugation tst''^ from {e, r^cr^, r^^crj;j,,T^crj.g2} to {e, r^zCr^, T^zCfxy, (^xyz} = 
{Txz<^z,crxyz) which is the Rxz symmetry group of GHC09. Of these symmetries, the 
^-localization breaks TxzO'z and TxzC^xy, since they involve periodicity in z, and leaves 
only axyz symmetry, which is in fact the symmetry of the localized equilibrium of plane 
Couette reported in SGBIO. For choices of z phase that are not integer multiples of iz/4:, 
each of the three symmetries of the periodic solution is broken by the localization, leaving 
completely unsymmetric initial guesses, corresponding to the rung solutions of SGBIO. 

2.5. Spanwise localized equilibria of plane Couette flow: computation 

In this section we construct new spanwise localized equilibrium solutions of plane Cou- 
ette flow by applying windowing and refinement to the spatially periodic EQ7 solu- 
tion. Figure [2] illustrates how spanwise-localized initial guesses with different symmetry 
groups are constructed from different spatial phases of the spatially periodic solution. Fig- 
ure l2ja) shows EQ7 at Re = 400 with fundamental streamwise and spanwise wavenum- 
bers a,7 = 1,2 (i.e. periodic lengths £x = 27r/a = 2tt and iz = 27r/7 = tt). The figure 
shows three copies of the periodic structure in the z G [37r/2, 3tt/2] subset of the full 
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Lx, Lz = 27r, Gtt computational domain. The spatial phase of EQ7 in figure ^a) is chosen 
so that one concentrated vortex structure is centered on the z = plane. In this phase 
the solution has symmetry group {o'xy,Txz,TxCrz)- Each of these symmetries is readily 
apparent in the figure, keeping in mind that the orientation of swirling with respect to 
the X axis and thus color of the isosurfaces changes under axy (rotation about z axis) 
and Txfyz {x shift, z-reflect symmetry). The windowing function W{z) is plotted as a 
function of z in the front face of the box. Multiplication of the periodic structure shown 
in figure p[a) by the windowing function, followed by projection onto the divergence- free 
subspace, produces the initial guess shown in figure l2jc) . In this case the windowing 
parameters a, 6 = 0.3, 1 were chosen to preserve the single concentrated vortex structure 
centered on the z = plane and to taper rapidly to nearly zero before the next vortical 
structure. As discussed in § |2.3[ localization in z breaks the Txz symmetry of the periodic 
solution, since it involves z periodicity, leaving a localized initial guess with symmetry 
group {(Jxy^TxCTz)- Any solution in this symmetry group will be an equilibrium, since the 
X reflection in axy prevents travelings waves in x and the z reflection in TxOz symmetry 
prevents traveling waves in z. 

Figure ^b,d) illustrate construction of a localized initial guess with different sym- 
metries by windowing the periodic solution in a different spanwise phase. Figure [2|b) 
shows the same periodic EQ7 solution as in figure 0[a) , but translated by a quarter- 
wavelength in z. In this phase the periodic NBCW solution has symmetries {(Txy, <Jz, Txz) ■, 
which are again readily apparent in the figure (and which can be derived by conjugating 
{(^xy^TxzTTx(Jz) with Tz and choosing the specified symmetries as generators for the 
conjugated group). A wider windowing function, with a = 6 = 1, preserved a pair of 
mirror-symmetric concentrated vortical structures in the core region and tapered rapidly 
to nearly zero before reaching the next vortices, as shown in figure ^b). The windowing 
breaks symmetries with factors of Txz , producing an initial guess for a spanwise-localized 
solution with symmetry group {oxy^iJ:^- Again, sign changes in all three coordinates in 
this symmetry group fix the phase of the velocity field with respect to the origin and 
rule out traveling waves. Thus, the localized solutions in both choices of z phase will be 
equilibria, unlike the localized solutions of SGBIO, where one choice of phase produces 
equilibria and the other streamwise-traveling waves. 

The localized initial guesses depicted in figure l2[c,d) were then refined to numerically 
exact equilibrium solutions of plane Couette flow shown in figure l3Fa,b), using a Newton- 
Krylov-hookstep search algorithm. The search algorithm finds equilibria as solutions 
of the equation f"^(u) — u = 0, where f"^ is the time-T integration of the discretized 
Navier-Stokes equations with appropriate boundary conditions. The time integration is 



performed with a Fourier-Chebyshev-tau scheme in primitive variables (Spalart et al. 



(1991); Canuto et al. (2006)) and 3rd-order semi-implicit backwards differentiation time 
stepping (Peyret (2002)). The Fourier-Chebyshev spatial discretization of u(x,i) takes 
the form 

J K L 

j = -J k=-K 1=0 

Spatial discretization levels are specified by (J, K, L) or by the corresponding iV^, x iVy x 
Nz = (2J-)-l) X L X {2K+1) physical grid. Nonlinear terms are computed with 3/2- 
style dealiasing. We set spatial discretization levels so the maximum truncated Fourier 
and Chebyshev modes are O(10~^) and O(10^^°) respectively. Experience has taught 
us that coarser spatial discretizations can result in spurious solutions. Symmetries were 
enforced through the search by projecting u — > (u + ctu)/2 for each of the generators a 
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of the appropriate symmetry group at the intervals AT — 1 during time integration. The 
residual of the search equation is ||f"^(u) — u||/T using the L? norm 

-,1/2 



V jv 



u • udx 



(2.5) 



where V is the volume of the computational domain. We measure the accuracy of a given 
solution by increasing spatial resolution by a factor of 3/2 in each direction, decreasing 
the time step by a factor of 2, and then recomputing the residual. Further details of the 
implementation of the search algorithm and time integration are given in GHC09, and 



the code is available for download at www.chaiinelflow.org (Gibson (2013)). 

The localized plane Couette solutions depicted in figure [3][a,b) were computed at Re = 
400 in a 27r, 67r computational box and a 16 x 65 x 128 grid, with integration time T = 10 
and time step dt = 0.077, resulting in a CFL number of about 0.6. The residuals of 
the windowed initial guesses began at roughly 10"'^ and were reduced to 10"^^ after six 
or seven iterations of the Newton-Krylov-hookstep algorithm. The O(10~^) accuracy of 
these solutions, measured as discussed above, is the best one could expect given that 
spectral coefhcients were truncated at that level. The tails of the localized solutions drop 
to O(10~'*) at the z — ±37r edge of the computational domain. The computational cost is 
modest: about one CPU-hour for each solution, running serially on a desktop computer 
with a 3.3 GHz Intel i7-3960X processor. 

2.6. Spanwise localized equilibria of plane Couette flow: roll-streak structure 

Figure [3] shows six spanwise-localized equilibrium solutions of plane Couette flow with 
streamwise wavenumber a = 1 and Re = 400. Figure p[a-d) were obtained by the 



localization and search methods outlined in § 2.5 (a,b) as described in detail in § 2.5 
and (c,d) by increasing the core region of the window to fit three and four copies of 
the concentrated vortex structures respectively. Solutions (a-d) show that it is possible 
to obtain localized versions of EQ7 with 1,2,3, and 4 copies of the basic concentrated 
vortical structure shown in isolation in (a) , by choosing appropriate centers of symmetry 
and width of the core region of the windowing function. We will refer to the solutions 
(a,b,c,d) as EQ7-1,2,3,4 respectively. Odd-numbered solutions have {(yxy,Tx(yz) symmetry 
and even-numbered solutions {a-^y, a^). We note that each of the EQ7-1,2,3,4 solutions lies 
on a distinct solution branch. This is in distinct contrast to the localized versions of the 
NBCW solutions, for which all solutions with the same symmetry lie on a single solution 
branch, and additional copies of the fundamental structure are added in a continuous 



fashion under continuation in Reynolds number, via homoclinic snaking (Schneider et al. 



(2010)). The localized EQ7 solutions, in contrast, can each be continued in Reynolds 
number around a single saddle-node bifurcation, and solutions on the opposite branch 
have distinctly different physical structure. For example, figure Isje) and (f) show the 
opposite branches of (b) EQ7-2 and (d) EQ7-4, obtained by continuing these downward 
in Reynolds, around a saddle-node bifurcation, and back upwards to Re = 400. 

The spatial structure of EQ7-1 and EQ7-2 is illustrated in more detail in Figure [4] 
Figure Sa-e) show the cross-stream velocity [v,w]{yz) for EQ7-1 in five streamwise- 
normal cross-sections spaced evenly between x — — tt and a; = 0, which are the front face 
and the middle of the boxes in figure[3|^a,b). The blue isosurface of signed swirling strength 
in the front half of the box in figure Isfa) appears here as a concentrated counter-clockwise 
vortex that begins just below and to the left of the origin at a; = — tt in (a), increases in 
strength and moves upward and to the right in (b-d) , and ends above and to the right of 
the origin at x = in (e). By the TxCFz symmetry of EQ7-1, the equivalent quiver plots 
for a; = through x — ir would be the z-mirror images of (a-e), showing a concentrated 
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Figure 3. Spanwise- localized equilibrium solutions of plane Couette fl ovy constructed 
from windowing the spatially periodic EQ7 equilibrium of Gibson et al. (20091 at Re — 400, in 
the z € [— 37r/2, 37r/2] subset of the Lx,Lz = 27r, Ott computational domain. (a,b) are localized 
solutions produced from the initial guesses shown in figure[2[c,d). See text for discussion of (c-f). 
The plotting conventions are the same as in figure [2] but additionally isosurfaces of streamwise 
velocity at u = 0.18 are plotted in red to show the positions of high-speed streaks near the lower 
wall (symmetric high-speed streaks near the upper wall are not shown). 



clockwise vortex starting below and to the right of the origin and moving upwards and 
leftwards, and corresponding to the green isosurface of signed swirling strength in the 
back half of the box in figure [sl a). Likewise, the predominant features of EQ7-2 shown 
in figure pTf-j) are two z-symmetric counter-rotating vortices that begin near the lower 
wall at a; = — tt in (f), and rise upwards and closer together in (g-j). Due to EQ7-2's axy 
symmetry, cross-sections from a; = to tt would appear as the y-mirror images of (j-f), 
with a pair of nearby counter-rotating vortices in the lower half of the {y, z) plane rising 
and separating, as in the back half of figure Islb). The clear and distinct concentrations 
of circulating fluid in figure |4] constitute our main justification for speaking of "regions of 
concentrated vortical structure" and show that isosurface plots of signed swirling strength 
are not misleading but rather show precisely where such regions lie. 

The mean roll-streak structure of EQ7-1 and EQ7-2 is illustrated in Figure[5] The four 
counter-rotating vortices surrounding the origin in figure p[a) result from x-averaging 
the counter-clockwise vortex that slopes upwards and rightwards from x = — tt to in 
figure p[a) and figure E[a-e) with its clockwise rj;CTz -symmetric counterpart that slopes 
upwards and leftwards from a; = to tt. These four vortices create the pattern of al- 
ternating positive and negative streamwise streaks (relative to laminar flow) shown in 
(b) by advecting high-speed fluid (utot = il) fi^om the walls towards the interior in the 
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Figure 4. Cross-sections of spanwise-localized equilibrium solutions of plane Couette 

flo'w. (a-e) Quiver plots of [v,w]{y,z) for EQ7-1 at a; = {— vr, — 37r/4, — 7r/2, — 7r/4, 0} with y 
vertical and z horizontal, (f-j) EQ7-2 at same x values, a = 1 and Re = 400 for both. Compare 
to the same solutions shown with isosurfaces of signed swirling strength in figure [3|a) and (b); 
the cross-sections here are evenly spaced along the front half of those figures. The z G [— tt, tt] 
subset of the full z £ [—Stt, Stt] computational domain is shown. 




(c) 
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(d) 

Figure 5. Streamwise-averaged roll-streak structure of spanwise-localized equilib- 
ria of plane Couette flow, (a) Quiver plot of x-average [v,'w\{y,z) and (b) contour plot of 
x-average u{y, z) for EQ7-1 with y vertical and z horizontal. (c,d) Same for EQ7-2. Both solutions 
are shown at a = 1 and Re = 400. Contour lines are plotted at levels u = ±{0.03,0.09,0.15}, 
with negative values in dashed lines positive in solid. Quiver plots are autoscaled. The z G [— tt, tt] 
subset of the full z € [— Stt, Stt] computational domain is shown. 
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region near z = 0, and low-speed fluid (utot ~ 0) from interior towards the walls for 
larger z. Figure p[c,d) shows the corresponding mean roll-streak structure for EQ7-2; 
here the doubling of the basic concentrated vortex structure compared to EQ7-1, appar- 
ent in figure ^b), results in the eight counter-rotating mean vortices shown in (c) and 
an increased pattern of alternating streamwise streaks shown in (d) . 



3. Traveling wave solutions of channel flow 

In this section we extend the results of §[2] to channel-flow conditions. For the first 
set of solutions, we use a numerical continuation method similar in spirit to Waleffe's 
homotopy of the NBCW solution between plane Couette and "half-Poiseuille" flow with 



no-slip at the lower wall and free-slip conditions at the upper wall ( Waleffe ( 1998 2001 )). 
When extended by symmetry to full channel conditions, Waleffe's continuation produces 
a traveling-wave solution symmetric about the channel midplane, with two NBCW-like 
roll-streak structures, each positioned in the high-shear region near either wall, and mirror 
symmetric (ay) to each other across the midplane. 

In the present study, we continue the EQ7 solutions from plane Couette to full channel 
conditions, enforcing no-slip boundary conditions on both walls throughout. The contin- 
uation is done in two stages: first with fixed wall speed and increasing pressure gradient, 
then fixed pressure gradient and decreasing wall speed. As in § |2.1[ we decompose the 
total velocity field into a base flow and a deviation, Utot (x, t) — u(x, t) + U{y) e^, and the 
total pressure field into ptot — p(x, i) + x dP/dx^ where dP/dx is a parametric constant 
corresponding to the externally imposed mean pressure gradient. To make the decompo- 
sition unique, we specify that the fiuctuation pressure p is periodic (so that the spatial 
mean of Vp is zero), and that the laminar velocity profile satisfies the no-slip conditions at 
the walls and balances the imposed mean pressure gradient, dP/dx = uU" , where v is the 
kinematic viscosity of the fiuid. Consequently the base flow is the laminar solution for the 
given viscosity, mean pressure gradient, and wall speed, and the fluctuation velocity satis- 
fies Dirichlet conditions at the walls. The Navier-Stokes equations again take the form of 



(2.1 ). The Reynolds number Re = Uh/v is based on a velocity scale U appropriate to the 
flow as it transforms from plane Couette to pressure-driven channel conditions, namely, 
U is half the relative wall speed when continuing in pressure gradient and the centerline 
velocity of the laminar base flow when continuing in wall speed [U — \dP/dx\ h? /{2v)). 
Thus in nondimensional terms the continuation is first in mean pressure gradient dP/dx 
from to —2/Re with wall speeds fixed at U{±1) — ±1 (equivalently from U{y) — y 
to U{y) ~ 1 + y — y'^), and then continuation in wall speed from 1 to with mean 
pressure gradient held fixed at P^ = —2/Re (equivalently from U{y) — 1 + y ~ y'^ to 
U{y) — 1 — y"^). For channel flow conditions U{y) = 1 — y^ and dP/dx ^ 0, eqn. (2.1 ) and 



boundary conditions are invariant under any combination of x and z translations and 
reflections about the y = and z = midplanes; thus the symmetry group of channel flow 

is {ay,az,T{Ax, Az)), where ay : [w, w, w]{x, y, z) -^ [u, —v, 'w\{x, —y, z). 

3.1. Spatially periodic traveling wave solutions of channel flow 

TWl: Figure [6](a) and figure [t]^ a, d) show a spatially periodic travehng-wave solution of 
channel fiow with symmetry {(Ty,az,Txz) obtained by continuation from plane Couette 
conditions, as described above. The starting point for continuation was the spatially 
periodic EQ7 equilibrium of plane Couette fiow at Re — 2000 and a, 7 — 1,2; the TWl 
channel traveling wave has the same spatial periodicity and is shown at Re = 2300. Note 
that TWl is mirror symmetric about the y = midplane. The y mirror symmetry is most 
clearly seen in the streamwise-averaged plots of cross-stream and streamwise velocity in 
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Figure 6. Spatially periodic traveling waves of channel flow w^ith decreasing wal- 
l-normal symmetry, (a) TWl, constructed by continuing EQ7 from plane Couette to chan- 
nel conditions, has symmetry group {(Jy,a^,T^z). Isosurfaces of signed swirling strength are at 
s = ±0.04 are shown in green/blue. High-speed streaks near the lower wall are shown by iso- 
surfaces of streamwise velocity at u = 0.02. High-speed streaks near the upper wall, symmetric 
to those near the lower wall, are not shown as isosurfaces but they are indicated in the contour 
plots of streamwise velocity in the back plane, (b) TW2, found from an initial guess judiciously 
chosen from numerical simulation data, has symmetry group {az,Txz)- Isosurfaces of signed 
swirling strength at s = ±0.10 are shown in green/blue, and high-speed streaks near the lower 
wall are shown by isosurfaces of streamwise velocity at u = 0.03. There are similar but weaker 
vortex structures and high-speed streaks near the upper wall, but they do not appear at these 
levels for the isosurfaces. Both (a) and (b) are shown at 0,7 = 1,2 and Re = 2300. (c) TW3, 
constructed by continuation in Reynolds number from TW2 and further symmetrization and 
Newton-Krylov refinement, has {a^, Txz) but periodicity a, 7 = 1, 6, and is shown at Re = 4000. 
Isosurfaces of signed swirling strength are at s = ±0.9 and it = 0.03. Each back plane shows 
contour plot of streamwise velocity at levels u — {0.03, —0.03, —0.09, —0.15, ...}, positive in solid 
lines and negative in dashed. 




wmm 



SSy^SlDj^^^'' ^^^^^V'' 



M 




b'^^-^cZD^--- d 



- 


--- 


-- 





-- 





- 


1 


^ 1 


'": 




^ > 


'^: 


/ 


^ 






^ - 






, 


) 


'O 


-0 











■c 



(f) 



ir/2 



Figure 7. Streamw^ise-averaged roll-streak structure of spanwise-periodic traveling 
waves of channel flow. Quiver plots of a;-average [v,w]{y, z) and contour plots of x-average 
u(y, z) for (a,d) TWl, (b,e) TW2, and (c,f) TW3, with the same spatial and Reynolds parameters 
as in figurelm Contour lines are shown at levels u — —0.2, —0.1 (dashed) and u = 0.025 (solid), 
with an additional negative contour line is shown aX u — —0.3 in (e). The quiver plots are 
autoscaled. Axes are z horizontal and y vertical. 



figure P7ia,d). Although ay symmetry is within the symmetry group of channel flow, we 
did not expect the continuation to produce a solution with this symmetry, since it was 
neither present in the initial plane Couette solutions nor allowed in the intermediate 
steps in the continuation from plane Couette to channel conditions. Instead, we expected 
that the increasing y asymmetry under continuation in pressure gradient would push 
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the vortex structures towards the lower wall, where the shear of the base flow is higher 
(|C/'(— 1)1 = 3, compared to |C^'(1)| = 1, for the base flow U{y) = l+y—y^ attained at the 
end of the pressure continuation) , and that this y-asynimetry would be maintained during 
continuation in wall speed down to U{y) = 1 — y^. However, it turned out that weak 
vortices formed near the upper wall under pressure continuation and grew in strength 
during wall-speed continuation so that the solution gained ay symmetry when the wall 
speed reached zero. In terms of the symmetry groups, the starting plane Couette solution 
had symmetry group {o'xy,az,Txz), continuation to nonzero dP/dx broke axy symmetry, 
and ay was gained at the final step of wall-speed continuation, resulting in symmetry 
{ay,az,T,j;z). The structure of vortices and streaks in TWl can roughly be described as 
two copies of EQ7 stacked on top of each other, with the upper copy either phase-shifted 
by half a wavelength in x or having [v, w] reversed in sign via mirror symmetry in y. This 
is apparent from comparison of TWl in figure l6t a) to EQ7 in figure llTb) and EQ7-1 in 
Hb). 

TW2: Figurel6[b) and figurel7Fb,e) show a spatially periodic traveling-wave solution of 
channel flow with asymmetry in y and symmetry {a^, Txz)-, obtained from an initial guess 
from turbulent simulation data ([Viswanath^ (2007]) ; Gibson et a/.^(2009J). Specifically, we 



z-mirror-symmetrized an arbitrary turbulent velocity field of channel flow at Re ~ 3750 
in a 27r, Itt box by applying u — >■ 1/2 (1 -I- CTz)u and then quenched the turbulent field 
by lowering the Reynolds number and continuing time integration with the bulk velocity 
fixed at 2/3 and the az symmetry enforced by projection at regular intervals. After some 
experimentation, we found that after quenching to Re = 2650, the fine-scale structure of 
the velocity field and the spatial-mean wall shear decreased quickly, the latter reaching 
a local minimum after about 50 time units and growing slowly again for another 50 time 
units before resuming a high level of wall shear with rapid fluctuations. The smoothness 
and length of this minimum suggested a close pass to a hyperbolic edge state. Using a 
velocity field from this minimum as an initial guess for a Newton-Krylov search produced 
a numerically exact spatially periodic, wall- localized traveling wave solution of channel 
flow with {azTTxz) symmetry. 

The y asymmetry of TW2 is exaggerated in figure ^Vo) by the binary character of 
isosurface plots. In fact at Re — 2300 TW2 has weaker vortex structures near the upper 
wall with swirling strength comparable to those of TWl, and weaker streaks there as 
well. Both the vortices and streaks near the upper wall are visible in the streamwise- 
average plots of TW2 in figure [7lb,e). Note also that the swirling strength isosurfaces 
of TW2 in figure [6lb) are at s = ±0.10, over twice the magnitude of those for TWl in 
figure IBFa) at s = ±0.4. Compared to the four layers of counter-rotating mean vortices 
stacked symmetrically about y = in TWl (see Figure l7[a,d)), TW2 has two layers of 
counter-rotating mean vortices below y = and a single layer of vortices above y = 0, 
and these vortices have the same orientation as the vortices below them. And though it is 
not clear from the autoscaled quiver plots, the mean vortices of TW2 near the lower wall 
are about three times the magnitude of those of TWl, as measured by magnitudes of the 
[w, w] velocities, and the mean vortices of TW2 above y = are of comparable magnitude 
to those of TWl. However, the y-asymmetry of TW2 increases as the Reynolds number 
is increased (see TW3) and as z-periodicity is relaxed (see TW2-1 and TW2-2). 

TW3: Figure l6Fc) and figure [7tc,f) show a wall-localized, spanwise and streamwise 
periodic traveling wave of channelflow with symmetry (ctz, t^z), discovered through con- 
tinuation of TW2 in Reynolds number. The fundamental z wavenumber of TW2 is 7 = 2, 
but as as Re increased towards 4000, the structure at this wavenumber weakened and 
structure at 7 = 6 grew, while the structure away from the lower wall weakened substan- 
tially until it became nearly laminar. TW3 was computed by zeroing all modes in TW2 
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Figure 8. Spanwise-localized traveling wave solutions of channel flow, (a) TWl-1, 
(b) TWl-2, (c) TW2-1, and (d) TW2-2. (a,b) TWl-1 and TWl-2 are {ay.T^a^,) and {ay, a,) 
symmetric traveling waves obtained by localizing TWl in two different z phases. (c,d) TW2-1 
and TW2-2 are r^a^ and CTz symmetric traveling waves obtained by localizing TW2 in two 
different z phases. Plotting conventions are the same as in figurelGJ but with isosurfaces of signed 
swirling strength and streamwise velocity at (a,b) s = ±0.05, u = 0.05, and (c,d) s — ±0.09, 
u = 0.08. Contour lines of streamwise velocity are shown in the back y, z plane at levels u — 0.03 
(solid lines) and (a,b) u = {-0.03, -0.09, -0.15}, (c,d) u = {-0.03, -0.09, -0.15, -0.21, -0.30} 
(dashed). In (b) a high-speed streak near the upper wall, symmetric to that near the lower wall, 
is suppressed to avoid visual clutter. By contrast, (c,d) show true asymmetry in high-speed 
speed streaks: (c) has high-speed streaks near the lower wall only, and in (d) the streaks near 
the upper wall are substantially weaker, below the given isosurface levels. Solutions are shown 
at Re = 2300 in 2: G [— 37r/2, 37r/2] subsets of their full Lx,Lz = 27r, 6n computational domains. 



at Re = 4000 with 7 < 6 and refining this initial guess to an exact traveling wave with 
Newton-Krylov-hookstep search. The resulting TW3 solution has periodicity a, 7 = 1,6 
and symmetry {(Jz,Txz), where t^z is understood as involving a half-cell shift in z with 
respect to the smaller iz = tt/S periodic length. Its most notable property is its very 
strong localization in the wall-normal direction, as evidenced by figure l7lc,f). 

3.2. Spanwise localized traveling wave solutions of channel flow 

In this section we construct spanwise-localized traveling wave solutions of channel flow 
by windowing the spanwise-periodic travelings waves of § |3.1| TWl and TW2 in different 
spatial phases. The resulting solutions are illustrated in figure [8) figure |9] and figure [TO] 
TWl-1, shown in figure IsFa), was formed by phase-shifting TWl in z by 7r/4 to 
give it symmetry group {<yy,(7zTx,Tj.z), extending this periodic solution from a 27r, Itt 
to a 27r,67r box, windowing the extended periodic solution, and then applying Newton- 
Krylov-hookstep refinement to the windowed initial guess, producing a spanwise-localized 
traveling wave with (cj,, CTzT^;,) symmetry. The windowing parameters a,b = 0.5,0.4 
were chosen to isolate the central vortex structures in swirling-strength plots. It did not 
take a great deal of effort to find windowing parameters that gave a successful initial 
guess; we merely adjusted a and b by tenths until swirling-strength plots of the initial 
guess appeared similar to the resultant solution in figure [sla). The initial guess for 
the wave speed was set to the wave speed of the underlying periodic solution TWl, 
and the pressure gradient was held fixed at dP/dx = —2/Re where Re = 2000. The 
Newton-Krylov-hookstep search converged to machine precision in six steps, consuming 
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Figure 9. Cross-sections of spanwise-localized, near-wall traveling waves of channel 

flow, (a)-(e) Quiver plots of [«,«;] (y, 2) for TW2-1 at a: = {-tt, -O.Gtt, -0.27r, 0.27r, 0.67r}, (f)-(j) 
same for TW2-2, with y vertical and z horizontal. Both solutions are at q = 1 and Re = 2300. 
The z G [— TT, n] subset of the full z G [— Stt, Stt] computational domain is shown. 
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Figure 10. Streamwise-averaged roll-streak structure of spanw^ise-localized traveling 
waves of channel flow, (a) Quiver plot of a;-average [v,w] and (b) contour plot of a;-average 
u{y,z) for TW2-1. (c,d) Same for TW2-2. a = l,Re ^ 400 for both solutions. Contour lines 
are plotted at levels u — {—0.25, —0.15, —0.05, 0.05, 0.15}, with negative values in dashed lines 
and positive in solid. Quiver plots are autoscaled. The z G [— tt, it] subset of the full [— Stt, Stt] 
computational domain is shown. 
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about half an hour of single-core CPU time on the machine described in § |2.5[ for a 
20x81x256 discretization of the 27r,67r computational domain. The computed solution 
turned out to be smoother in x and z than the periodic solution on which it was based, 
enough that 10~^ x 10^^" x 10"^ truncation levels were retained on a reduced grid of 
12x81x192. Similarly TWl-2 figure IsFb) was found by refining an initial guess formed 
from windowing TWl in in its z-phase with (ctj,, az, Txz) symmetry and with windowing 
parameters a,b — 1.2, 0.4, resulting in a traveling wave solution with {ay, a^) symmetry. 
TW2-1 and TW2-2, shown in figure IMc,d), are likely the most interesting solutions 
presented in this paper, as they represent traveling wave solutions that are spanwise 
localized and strongly concentrated near a single wall, and as such are the solutions 
most likely to be relevant to the lambda-vortex coherent structures that occur near the 



walls of shear flows (Saiki et al. (1993)). TW2-1 was formed by windowing TW2 in its 
z-phase with symmetry {txUziTxz) and windowing parameters a,b = 0.6,0.4 to get an 
initial guess with r^az symmetry, and refining that with Newton-Krylov-hookstep. We 
were unable to form TW2-2 by the phase-shifting, windowing, and refining procedure 
employed for other localized solutions. Instead TW2-2 was formed by shifting TW2-1 
leftwards in z until its right-hand high-speed streak was centered on z = 0, extending the 
shifted field on z ^ onto z ^ by z-mirror symmetry, and then refining this initial guess 
with Newton-Krylov-hookstep. We note that TW2-1 and TW2-2 required finer grids to 
resolve the stronger vortex structure adequately compared to TWl-1 and TWl-2. Larger 
computational domains were also necessary, due to the increased width of the streamwise 
velocity deficit (visible as contour lines in the back plane of figure Isl) . The computations 
of TW2-1 and TW2-2 from initial guess to took about seven CPU-hours each. 

Each of the traveling wave solutions depicted in figure [8] is formed from variations of 
very similar basic structure, which mostly clearly seen in isolation in TW2-1 in figurepTc). 
Near the lower wall there is an x-periodic chain of concentrated vortices, alternating in 
sign of circulation (clockwise/counter-clockwise) , and nearly aligned with the x-axis but 
tilting slightly in the wall- normal and spanwise directions, as shown in figure ^a-e). 
The tilting of the chain of alternating vortices results in a nonuniform x-average in the 
cross-stream fiow, specifically, a pair of counter-rotating mean vortices near the wall, 



figure 10 ^a). The mean vortices draw low-speed fluid upwards between them and high- 
speed fluid downward on either side, producing the mean high-speed streaks on either side 
of the mean vortex pair, depicted by solid contour lines in figure lIoFb). The streamwise 
momentum exchanged induced by the near-wall mean vortices has a net negative effect: 
the region of mean streamwise flow slower than laminar indicated by negative (dashed) 



contour lines in figure 10 b) is larger in both y, z area and magnitude than the high-speed 
streaks that outpace laminar flow, indicated by solid contour lines. Thus the net effect of 
the roll-streak structure is a decrease of bulk flow relative to laminar for a fixed pressure 
gradient. TW2-2 roughly consists of two copies of TW2-l's basic structure, repeated with 
mirror-symmetry about the z = plane. This is evident from comparison of TW2-1 and 
TW2-2 in figuregc) and (d), figme^a-e) and (f-g), and figure [lO]; a, b) and (c,d). TWl-1 
and TWl-2 have similar structure to TW2-1 and TW2-0, but with mirror symmetry in 
y and weaker magnitudes of vorticity (by roughly a factor of two) and net velocity deficit 
relative to laminar (by a factor of four or more) . 

The cross-stream quiver plots of these exact wall-bounded traveling waves, figures |8] 
and [To] bear intriguing resemblances to the cross-sections of near- wall lambda vortices of 



developing channel flow shown in Saiki et al. (1993) figures 10 and 11. For example, the 



streamwise mean v,w fiow of TW2-1 in figure 10 is strikingly like the temporal-mean 



V, w of the K-type-generated lambda vortex shown in|Saiki et al. ( 1993 ) figure 10(b). The 



A-structure shown in Saiki et al. (1993) figure 2(b) is evident in TW2-2 in figurc^d). Of 
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Table 1. Characteristics of traveling- wave solutions of channel flow. 



course the traveling waves in this paper are exact solutions of different flow conditions 



than the spatially-developing channel flow of Saiki et al. ( 1993| , so some translation of 
the flow conditions is required to make more precise comparisons. We intend to pursue 
this in future research. 

Table [l] summarizes physical characteristics and discretization properties of the com- 
puted traveling-wave solutions of channel flow. Wave speed c is in nondimensionalized 
units where the center line velocity of laminar flow is U — 1. Lx,Lz are the size of the 
computational domain. Accuracy is measured as |jf^(u)) — ru||/r for T = 10, where 
f "^ is the time-T forward time integration of the Navier-Stokes as performed by direct 
numerical simulation, and where t = t{cT,0) is the x-translation resulting from the 
computed wave speed c over time T. The magnitude of the tails is given as the maximum 
of |w|,|i;| and \w\ over x,y at z = ±Lz/2, the edge of the spanwise-extended periodic 
box. 



4. Discussion 

4.1. Exponential decay of tails 

The spanwise localized solutions presented in this paper display a three-part structure: a 
core region that closely resembles a periodic solution, a transition region, and weak tails 
that decay to laminar flow. In this section we show that the tails of spanwise-localized 
streamwise-periodic equilibria are dominated by a mode that decays exponentially at 
e""!"^! , where a is the fundamental streamwise wavenumber, and that the structure of the 
tails depends on flow parameters a. Re, the laminar flow proflle U{y), and the wavespeed 
c, but not on the details of the solution's core region. 

As the tails of a spanwise-localized equilibrium or traveling wave approach laminar 
flow, the perturbation velocity u approaches zero, so we expect u to approximately 



satisfy the linearized form of (2.1 ) in which u • Vu is set to zero 



00 -1 

^ + f/^ + vf/'e, = -Vp+— V^u, V-u = 0. (4.1) 

We look for normal-mode solutions of the form 

u,- ^(x) = u,- ^(y) e'(^"(---*)+^-), p,- ^(x) = p,, ^(y) e^0"(---*)+7-) (4.2) 

where a is real and 7 = 7^ -I- iji has 7^ > for tails that decay exponentially as 2: — > 00. 
The slowest decaying normal mode solution, with the smallest positive %, will dominate 
the tails as z — )• 00. For the remainder of this section we drop the j,j subscripts. To 
eliminate pressure we convert to velocity-vorticity form by taking the y-components of 
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(b) 
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Figure 1 1 . Exponential decay in the tails of spanwise localized plane Couette equi- 
libria, (a) Componentwise decay rates for EQ7-1 at a = 1/2, Re = 600: u,v,'w (solid lines) 
scale as e"°^, e~^"^, e""^ (dotted lines), (b) Decay of ||«;|| 
to e~"^ (dotted lines) at Re = 600 and several values of a. 



for EQ7-1 (solid lines) compared 



the curl and the curl of the curl of (4.1 ) 

d 



at + ^l;'^ 



Re 



d_ 
Ft 



U 



d_ 
dx 



Re 



(4.3) 



where rj — u^ — Wx is the wall-normal vorticity. Boundary conditions are r](x, ±1, z) 
and v{x, ±1, z) — v'{x, ±1, z) = 0. 







Eqn. 4.3 and boundary conditions permit a number of types of solution. Where sym- 
metries allow, the solution that dominates behavior in the tails of localized plane Couette 
equilibria and traveling waves turns out to be the trivial solution v = rj = 0. This so- 



lution in conjunction with the divergence-free condition and the ansatz (4.2) requires 
that Uxx = —Uzz and Wxx = —Wzz, which is satisfied by 7 = zLija and w = -iiijii, the 
different signs governing exponential decay in the different limits z — > ±c». The j = 1 
mode with 7 = ±ia and w = ±iu thus gives the slowest exponential decay rate {j = is 
ruled out since it does not decay and thus cannot be part of a z-localized solution) . The 



y component of (4.1 ) gives that dp/dy = so that p{y) = p is a, complex constant whose 
magnitude and phase are set by the pressure conditions at the edge of the transition 



region at some fixed value of z. The x component of (4.1 ) gives 



u" ~ iaRe{U 



c)u — iaRe p. 



(4.4) 



The boundary conditions u(±l) — set the homogeneous solution to zero, so that u 
and thus w are determined by the fixed value of p, the base flow profile U{y), and the 
parameters a, c and Re. They are independent of the structure of the core-region solution 
(except for differences of x phase between the ±z tails resulting from the symmetries 
of the solution) . Thus the v = r] = solution to the linearized Navier-Stokes equations 
contributes to the tails an exponentially decaying mode of form [u{y), 0, i{t(y)]e'"^^~'^*^~"^ 
as z — > 00. 

Figure [TT] confirms that over a wide range of a, the tails of EQ7-1 are dominated by 



u, w components that decay as e~"^. Figure 11 a.) shows the inf-norm of u, v and w as a 
function of z for EQ7-1 at a = 1. In this context ||u||oo(z) is the maximum of \u(x, y, z)\ 
over x,y as a function of z. As argued above, the magnitudes of the u and w components 
are equal in the tails and scale as e~"^. The higher-order e~^"^ scaling for v results from 
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>. oa^^-'-z--: 



(a) 




(b) 

Figure 12. (a,b) Asymptotic form of the tails of spanwise- localized equilibria of 
plane Couette flow. Contours of streamwise velocity u{x,y) at fixed z: (a) a sfice of EQ7-1 
for a = 1 and Re = 400, at z = 11, where ||u||oo ~ 10~®, and (b) the asymptotic v — ri = 
normal mode solution. Contours are plotted at ±[0.15, 0.45, 0.75] times the maximum of u, with 
negative values in dashed lines, (c) Exponential decay rate for three types of normal 
modes in the tails of plane Couette equilibria at Re — 400. The minimal 7^ as a function of 
Q i s plo tted for (D) v = r] = solutions, (o) solutions of (4.51 with v{y) = 0, and (>) solutions 
of (|4^. 



the quadratic nonlinear term in ui,y,'Wi,y that has been suppressed on the right-hand 
side of (4.11 for the ^0,7 equation in for j = 0. For figure flljb), we continued EQ7-1 
paramctrically in a and observed that the e^"^ scaling holds over the explored range of 



1/4 ^ a ^ 2. For further confirmation of the dominance of u = 77 = modes, figure |12[a) 
compares an x, y slice of streamwise velocity u(a;, y, z) from EQ7-1 at a fixed z to (b) the 
same as computed from (4.4 1. 



The linearized equations for the tails have solutions other than v = rj = Q\ to show 
these decay more rapidly than the v = -q = Q modes (for plane Couette flow), we reduce 
eqn.|4.3|to the eigenvalue equations 



ia{U — c)fj + ijU'v 
ia{U - c){i" - (a^ + -/^)ij) ~ iaU"v 



1 

Re 
1 

Re 



(ry"-(«2 + ^2)~), 



(4.5) 



({}"" - 2(^2 + ^^)y" + (a2 + ^2)2-)^ (4 , 



The latter is the time-independent form of the familiar Orr-Sommerfeld equation for 



three-dimensional disturbances. Equation (4.61 is independent of fj and the eigenvalues 
7 can be found numerically for given a, c. Re, and U{y). The i^U'v coupling term in 



(4.5) acts as a nonhomogeneous forcing, requiring particular solutions for fj{y) to match 
the eigenmodes of (4.6). Eigenvalues distinct from those found for (|4.6|) can be found 



by solving (4.5) with v{'y) = 0. Figure 12 |c) shows the minimal 7^ allowed by these two 
equations as a function of a for plane Couette equilibria at Re = 400. Note that w = 77 = 
modes have the smallest 7j for a in the range shown (a < 2); thus these modes dominate 
the behavior of the tails in all solutions with streamwise wavelength greater than tt. It 
should be noted, however, that streamwise constant (a — 0) modes can exist in domains 
of any length Lx , and thus might play a dominant role in short (large a) domains if the 
symmetries permit them. 



4.2. Asymptotic scaling of streamwise Fourier harmonics 

In this section we provide a numerical account of the large-Reynolds behavior of the peri- 
odic and localized solutions developed in previous sections. In particular we measure the 
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scaling of various streamwise Fourier components of the solutions with i?e, and we show 
the development of critical layers at large Re. These features are key to the asymptotic 



analysis of NBCW suggested by Wang et al. (2007) and developed into a complete the 



ory by Hall & Sherwin (2010). The main results are as follows. The streamwise Fourier 



components of EQ7 and solutions related to it, localized and in channel conditions, obey 
scaling laws similar to those of NBCW, albeit with different exponents and substantially 



different magnitudes, suggesting that Hall & Sherwin (2010 )'s asymptotic analysis could 



be carried over to the new solutions. EQ8 and the j/-asymmetric channel flow solutions, 
in contrast, do not fit the asymptotic scaling framework so cleanly. All solutions appear 
to have well-defined critical layers, however, and the critical layer is particularly simple 
for EQ7 and its localized counterparts. 



As suggested by Wang et al. (2007) and developed into a complete theory by Hall & 



Sherwin ( 2010 ), a reduced PDE system can be developed for the spatially periodic NBCW 



solution from an asymptotic analysis of its streamwise Fourier modes and the critical layer 



that develops at large Reynolds numbers. Wang et al. (2007) showed numerically that 



NBCW has a simplified, quasi-2D structure in the limit of large Reynolds numbers, with 
a balance between an 0(1) streamwise-constant streaks, 0{Re^^) streamwise-constant 
rolls, and an 0(i?e~'^'^) mode in the first (fundamental) stre amwise Fourier harmon ic, 

-1/3). [Hall fc SherwinI ([2OI0I) in 



which concentrates in a critical layer of thickness 0{Re~ 
turn developed an asymptotic theory for NBCW based on vortex-wave interaction that 
provides insight to the physics of how these components of NBCW balance, predicts their 
scaling exponents, and which reduces the computation of the solution from a 3D Navier- 
Stokes problem at large- i?e to a 2D PDE at Re = 1 coupled with a linear wave evolution 
equation. Specifically, the interactions of very small fundamental-mode streamwise waves 
within the critical layer generate nonzero mean stresses that cause jumps in the pressure 
and the normal derivative of roll velocity across the critical layer. The jump in roll shear 



drives the mean rolls, which in turn drive the mean streaks. Hall & Sherwin (2010) show 
that these effects balance to leading order in Re~^, and that the asymptotic theory also 
reduces computation of the 3D steady state at high Reynolds number to a simpler 2D 
calculation at unit Reynolds number coupled with a linear wave evolution equation. 



This reduced quasi-2D PDE model of Hall & Sherwin ( 2010 ) is of particular interest to 



us since a theoretical analysis of spanwise localization in solutions of Navier-Stokes should 
be easier to develop in the context of a reduced model. There is strong numerical evidence 
that a theory of localization in solutions of Navier-Stokes might be developed. |Schneider| 
et al. (2009) noted a remarkable resemblance between the x, y-averaged energy of the 



localized NBCW solutions and localized solutions of the Id Swift-Hohenberg equation 
found by Burke & Knobloch (2007). The similarity was made more remarkable by Schnei- 



der et al. 



(2010)'s demonstration that the localized NBCW solutions undergo homoclinic 



snaking under continuation in Reynolds number, just as the localized Swift-Hohenberg 
solutions do under continuations in their bifurcation parameter. For Swift-Hohenberg, 
homoclinic snaking of localized solutions is quite well understood theoretically via "spa- 
tial dynamics" (Burke & Knobloch (2007)). Time independence reduces the 4th-order 



Swift-Hohenberg PDE for u{x,t) to 4th-ordcr ordinary differential equation (ODE) on 
u{x), which can then considered as a 4-dimensional dynamical system where x plays the 
role of time. In this view, spatially periodic solutions correspond to periodic orbits of 
the spatial dynamics, and spatially localized solutions correspond to homoclinic orbits 
that start at the origin u = at i — > — (X), grow away from it along unstable direction, 
wander at finite amplitude for some time, and then reapproach u = at i — ^ cx) along 
a stable direction. Localized solutions display approximately periodic form in their core 
regions when the finite-amplitude excursion away from u — makes a number of cir- 
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Table 2. Scaling exponents for streamwise Fourier harmonics of equilibrium and 
traveling wave solutions of plane Couette and channel flow. Components of solutions 
u scale as || ■ \\oo ~ 0{Re~'^) for the given values of /i. 
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cuits in the neighborhood of an unstable periodic orbit of the spatial dynamics. Thus 
the close correspendence between localized Navier-Stokes solutions and localized Swift- 
Hohenberg solutions points to the possibility of a theoretical explanation of localization 
in invariant solutions of Navier-Stokes. It also points to the importance of a reduced 
PDE system describing the localized solutions, ideally to a Id system in the spanwise 
coordinate, to which the idea of spatial dynamics might be applied. We report on the 
scaling of streamwise Fourier components and the development of a critical layer in the 
plane Couette and channel solution because they are essential ingredients for developing 
such a reduced-order PDE model. 

The perturbation velocity field of equilibrium and traveling wave solutions can be 
expressed as a sum of streamwise (x) Fourier modes in the form 



u(x,i) 



i 



<^jiy^z] 



,'ye 



(4.7) 



where = a{x — ct) and where u_j(j/, z) = u*(y, z) so that u is real- valued. The stream- 
wise constant mode Uq can be decomposed into streamwise streaks UQ{y,z) and cross- 
stream rolls [0, Vo, WoKy, z). Recall that u is the deviation from laminar flow U(ii)e^, with 
Utot = u + U{y)ex and Utot.o = uo + U{y), so that the streaks are defined relative to lami- 
nar flow. Wang et al.\\2007 ) define streaks relative to the z-averaged mean flow, but we do 
not, since z-averaging is inappropriate for spanwise-localized solutions. |Hall fc Sherwin 
( 2010 ) refer to the first (fundamental) harmonic Ui {y, z) or Ui {y, z) exp{ia{x — ct)) + c.c. 
as the wave mode. 

Figure [13] shows the scaling with Reynolds number of the magnitude of the streaks, 
rolls, waves, and second and third harmonics for a number of plane Couette and channel 
solutions. The magnitude is computed with the inf-norm; for example, the magnitude 
||ui||oo of the wave is the maximum over x^y,z and the vector components [u,w,w] of 
u = Ui(2/, z)e'"^ + u_i(2/, z)e~'"^. A number of these solutions exhibit very clear scaling 
^ OiRe- 



of the form llu,- 



-Pj 



For example Figure 13 



shows that for NBCW, the 
and the waves T7(i?e~°'^^). The streak and roll 



streaks are 0(1), the rolls 0{Re ^). 

values scalings equal the theoretical predictions of Hall & Sherwin (2010), and the wave 



scaling is within 2%. Note that the choice of inf-norm changes the scaling exponent for 
the NBCW wave component compared to the value of ^i = 0.9 reported by | Wang et al.\ 
pOOTl). 



EQ7 in figure 13 h) shows equally clear asymptotic scaling with exponents comparable 



but not equal to those of NBCW. The same is true of all solutions derived from EQ7 by 
parametric continuation and localization by windowing. Examples of such EQ7-related 



solutions are shown in figure 13 namely (d) TWl, the spatially periodic traveling wave 



24 



J.F. Gibson and E. Brand 



lo" 




(a) 



10' 



10 10 

Re 



(d) 



(g) 



(b) 




(e) 



(h) 




Figure 13. Scaling of strearti'wise Fourier modes of plane Couette and channel so- 
lutions. (a,b,c) Spatially periodic plane Couette equilibria NBCW, EQ7, EQ8, at q:,7 = 1,2. 
(d,e,f) Spatially periodic channel flow traveling waves TWl, TW2, TW3, at a,7 = 1,2. (g,li,i) 
Spanwise localized solutions: EQ7-1 plane Couette equilibrium and TWl-1 and TW2-1 channel 
flow traveling wave, at a = 1. The magnitude of various Fourier components of the velocity, 
as measured by inf-norm is plotted against Reynolds number. The labels s,r,l,2,3 indicate the 
streaks uo, rolls [i)o,wo], and the 1st, 2nd, and 3rd streamwise Fourier harmonics ui, U2, and 
U3. The vertical axis is the 00-norm of the given component (see text). 



of channel flow obtained from EQ7 by continuation; (g) EQ7-1, the spanwise localized 
equilibrium of plane Couette flow obtained by windowing; and (i) TW-1, the spanwise 
localized traveling wave of channel flow obtained by windowing TW-1. Scaling exponents 
for these solutions plus EQ7-2 and TW7-2 are listed in table[2J It should be noted that the 
magnitudes of Fourier harmonics of EQ7 and solutions derived fronr it are substantially 



different from those of NBCW. For example, comparison of EQ7 in figure 13 b) to NBCW 
in (a) shows that EQ7's streaks are about a factor of three smaller than NBCW's, and its 
fundamental harmonic is about a factor of three larger, resulting in an order of magnitude 
less scale separation between these components. 

In contrast, the upper-branch and y- asymmetric solutions reported here have no clear 
asymptotic scaling in streamwise Fourier harmonics and substantially poorer separation 
of scales, namely figure [l3J[b) EQ8, the upper branch of EQ7; (e) TW2, a spanwise 
periodic traveling wave of channel flow obtained from judiciously chosen DNS data; (f) 
TW3, a higher- wavenumber spanwise periodic traveling wave of channel flow obtained 
from continuation in Re of TW2; and (i) TW2-1, a spanwise localized traveling wave 
of channel flow obtained from windowing TW2. TW2-2 (not shown) is similar to TW- 
1. Among these, none of EQ8, TW2, TW2-1, or TW-2 continue in a straightforward 
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Figure 14. Critical layers of spanwise periodic equilibria of plane Couette flow. 

(a,b) NBCW at Re = 1000 and 30000, (c,d) EQ7 at Re = 1000 and 30000, and (e,f) EQ8 at 
Re = 1000 and 3000 for streamwise, spanwise wavenumbers 0,7 = 1,2. Dashed contour lines 
show total streamwise velocity at levels Utot ~ ±{0,0.25,0.50,0.75}. The critical layer where 
Wtot (y, «) = c = is shown with a thick solid contour line. Thin solid contour lines show the 
autoscaled magnitude of the fundamental Fourier harmonic, |ui|. The horizontal and vertical 
coordinates are 2, y respectively. 

fashion to higher Reynolds number; instead the solutions curves turn around at finite Re 
and follow complex paths. The same is true for EQ2, the upper branch of NBCW. The 
numerical evidence thus weighs against the possibility of an asymptotic analysis of these 
solutions based on streamwise Fourier harmonics. 



4.3. Critical layers 
The development of critical layers is an important consequence of the separation of scales 



in the streamwise Fourier modes (Wang et al. (2007) and Hall & Sherwin (2010)). The 



critical layer is the surface on which the mean streamwise fluid velocity matches the 
wavespeed, i.e. Mtot,o(2/, -z) = c. When higher harmonics become negligible and the roll 
velocities vq,wq are small compared to the streaky streamwise velocity '«tot,Oi the equation 
for the fundamental mode simplifies to 



[ia{ut. 



,ot,0 



c)ui + (ui • Vutot,o)e.] e'^ = V(pie'') + i?e-iv'(uie'*) 



(4.8) 



As argued by Wang et al. (20071, for large Re the fundamental harmon ic U i concentrates 
in a region of thickness 5 = Re~^^^ about the critical layer, in which (4.8) is dominated 
by a balance between its first and last terms. For a point x in this region and Xc nearby 
on the critical layer, this requires a balance of a{utot,o{^) — c) « a(x — x^) • Vutot.o 
against Re~^V^ . If (5 ^ |x — xd is the thickness of the region, the balance requires 
a^lVutot.ol ^ Re-^5-^ or 5 = (a|Vutot,o|i?e)-V^ 



Figure 14 illustrates the development of the critical layer for three spatially periodic 
equilibria of plane Couette flow. For equilibria, the wavespeed vanishes, so the critical 
layer in these plots is the surface y — f{z) on which utot{y,z) = 0. For NBCW, shown 
in figure fl4la,b), the height of the critical layer varies in z, and its thickness S decreases 
as Re~^/^etween (a) Re = 1000 and (b) Re = 30000. Figure ll4Fa,b) largely duphcates 



figures 2 and 3 of Wang et al. ( 2007 ) ; however we note that our plots show contourlines 
of |ui| rather than just the vertical component, \vi\, and so more clearly convey the fact 
that, for NBCW, the concentration of the fundamental mode is spread almost uniformly 
over the entire critical layer. 



26 



J.F. Gibson and E. Brand 



(b)-^ 



(a) -^-3 =2-10 1 2 3 (c) '^ 



(d)-^ 



'--''''''''^''(^''' 




^vH::^^;l^;;;;; 




-3-2-10 1 2 


3 
















Figure 15. Critical layers of spanwise localized equilibria of plane Couette flow. 

(a,b) EQ7-1 at Re = 1000 and 30000 and (c,d) EQ7-2 at Re = 1000 and 30000 for streamwise 
wavenumber q = 1, with y vertical and z horizontal. Plotting conventions are the same as 
figure 14 The z G [— tt, tt] subset of the full [— 37r, Stt] computational domain is shown. 
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Clr,y-'-3 -2-10 1 2 

Figure 16. Critical layers of spanwise localized traveling waves of channel flow. (a,b) 
TWl-1 at Re = 2000 and 30000 and (c,d) TW2-1 at Re = 2000 and 20000 for streamwise 
wavenumber a = 1, with y vertical and z horizontal. Plotting conventions are the same as fig- 
: except total streamwise velocity contours are shown at levels Utot = {0.1, 0.3, 0.5, 0.7, 0.9}. 
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The z e [— TT, it] subset of the full [— 37r, 37r] computational domain is shown. 



The critical layers for EQ7 and its upper branch EQ8 are markedly different (fig- 
ure fT4lc,d) and (e,f)). First, the critical layer for these solutions is the line y = 0. This 



follows from their Uxy symmetry, i.e. [u, w, w](x, y, z) = [— w,, — w, w](— a;, — y, z). Under 
this symmetry the x-average of the perturbation velocity and total velocity vanishes on 
y = 0. We note that some of the complexity of Hall & Sherwin (2010)'s analysis results 



from the need to work in a coordinate system aligned with the curved criticial layer; for 
EQ7 and EQ8 this complexity would be eliminated. Second, the fundamental mode Ui 
concentrates not uniformly over the whole critical layer, but apparently on isolated spots 
within it. Third, EQ8 seems to form a critical layer at Re — 2000, even though its scale 
separation is much poorer and its large- i?e limit apparently does not exist. 

Figure [15] shows that the critical layer structure of EQ7 carries over directly to its 
spanwise-localized counterparts EQ7-1 and EQ7-2, with tapering to laminar flow at large 
\z\. In particular the isolated concentrations Ui on the critical layer can be seen, by com- 
parison with figure l4l to result from the first-harmonic x variations of the y, z-localized 
and concentrated vortex structures. Figure [l6|^a,b) shows that EQ7 and EQ7-l's critical 
layer structure carries over to TWl-1, in two copies mirrored symmetrically about y = 0. 
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The channel travehng waves have a nonzero wavespeed c and lack EQ7's axy symmetry, 
and therefore have a critical layer wtot,o(2/j z) — c = whose height varies in z. 

Figure [l6[ c,d) shows critical layer development for the y-asymmetric channel traveling 
wave TW2-1. Note that y- asymmetry increases while y and z length scales decrease with 
increasing Re. In particular, Ui concentrates in a smaller region that approaches the 
wall as Re increases. An obvious question is whether this represents a near-wall coherent 
structure that is constant in wall units. We intend to address this question in future 
work. For the time being we note that the behavior illustrated in figure [T6J[c,d) is still 
subject to a prescribed length scale in the form of the streamwise wavelength a, and that 
this prescription must be removed, by streamwise localization or proper scaling with Re, 
in order for the lengthscales to be determined naturally. 

5. Conclusions 

We have found a number of new spanwise-localized equilibrium solutions of plane 
Couette flow and traveling- wave solutions of channel flow, and additionally a few spanwise 
periodic solutions of channel flow incidental to the construction of the localized solutions. 
The spanwise localized solutions consist of a core region that closely resembles a spanwise 
periodic solution, a transition region, and exponentially decaying tails. The decay rate 
of the tails is e^"'^', and their structure is determined by solely by the streamwise 
wavenumber, the laminar flow profile, and the wavespeed, and is otherwise independent 



of the structure of the core region. The solutions related to Itano & Generalis ( 2009 1 and 



Gibson et al. (2009 1 's HSV/EQ7 display clear scale separation and asymptotic scaling 
in streamwise Fourier harmonics, suggesting that they are amenable to analysis via a 
reduced-order PDE retaining only a few harmonics. 

Several solutions, namely TW2-1 and TW2-2, capture particularly isolated and el- 
emental exact coherent structures in the near-wall of shear flows, which suggestively 
resemble the lambda vortices that form in spatially developing flows. These solutions 
consist of long bands of concentrated vortices near the walls, with alternating orienta- 
tion, and roughly aligned with the streamwise axis but tilting slightly in the spanwise 
and wall- normal directions. The concentrated vortices near the walls are flanked by small 
high-speed streaks very near the walls and otherwise surrounded by very large regions 
where the streamwise velocity is reduced relative to the laminar background. These so- 
lutions capture, as exact time-independent solutions of Navier-Stokes, the process by 
which near-wall vortices exchange momentum between the near-wall and core regions of 
wall-bounded shear flows and thereby increase drag. 
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